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Electron spin resonance in a dilute magnon gas 
as a probe of magnon scattering resonances 
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We study the electron spin resonance in a dilute magnon gas that is realized in a ferromagnetic spin 
system at low temperature. A quantum cluster expansion is developed to show that the frequency 
shift of the single-magnon peak changes its sign and the linewidth reaches its maximum across a 
scattering resonance between magnons. Such characteristic behaviors are universal and can be used 
to experimentally locate the two-magnon resonance when an external parameter such as pressure 
is varied. Future achievement of the two-magnon resonance may have an impact comparable to 
the Feshbach resonance in ultracold atoms and will open up a rich variety of strongly correlated 
physics such as the recently proposed Efimov effect in quantum magnets. We also suggest that the 
emergence of an Efimov state of three magnons may be observed with the electron spin resonance 
through an additional peak structure in the linewidth caused by the three-magnon resonance. 
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I. INTRODUCTION 

Scattering resonances play an important role in 
physics. In the field of ultracold atoms, Feshbach res- 
onances induced with magnetic field are used to control 
the interaction between atoms.— This unparalleled tun- 
ability of the interaction has been essential to the recent 
remarkable advance of the field and led to experimen- 
tal realization of a rich variety of physics, such as the 
BCS-BEC crossover in Fermi gases^i^ and the Efimov ef- 
fect in Bose gases.'* These strongly correlated phenomena 
in the vicinity of the scattering resonance are universal, 
i.e., independent of microscopic details. Also, an atom 
loss peak caused by the three-atom or atom-dimer reso- 
nance has been used as a signature of the emergence of 
an Efimov trimer4i^ 

Turning to the field of condensed matter, the scattering 
resonance can be induced between collective excitations 
in ferromagnetic spin systems (magnons), for example, 
by tuning the exchange coupling.^ Because the ferromag- 
netic coupling is sensitive to the structure of ions,^ it is 
experimentally possible to vary the exchange coupling 
significantly with pressure^ and bring the system close 
to the two-magnon resonance. Future achievement of 
the two-magnon resonance may have an impact compa- 
rable to the Feshbach resonance in ultracold atoms and 
will open up a rich variety of strongly correlated physics 
such as the recently proposed Efimov effect in quantum 
magnets.— But, what is an experimental signature of the 
two-magnon resonance? 

In this article, we study the electron spin resonance 
(ESR) in a dilute magnon gas as a probe of magnon 
scattering resonances. ESR is a powerful experimental 
technique to investigate the dynamics of interacting spin 
systemS)2iiiS while theoretical calculations of its spectrum 
are in general challenging. Here we develop a quantum 
cluster (or virial) expansion which uses a small parameter 
called fugacity z = e~^^ <C 1 related to the diluteness of 
the system to perform a systematic expansion, where A is 



a single-particle excitation energy and /3 = l/{k-BT) is the 
inverse temperature We first demonstrate the quan- 
tum cluster expansion in Sec. |ll] by computing a single- 
particle spectral function in a dilute Bose gas, which 
serves as the universal formula for ESR at low temper- 
ature and can be used to extract the scattering length 
between magnons. Then in Sec. IIIIl we apply the same 
method to ESR in a dilute magnon gas and show that 
the frequency shift of the single-magnon peak changes its 
sign and the linewidth reaches its maximum across a scat- 
tering resonance between magnons. Such characteristic 
behaviors are universal and can be used to experimen- 
tally locate the two-magnon resonance when an external 
parameter such as pressure is varied. Finally, Sec. IIVI 
is devoted to the conclusion of this article and discus- 
sion where we also suggest that the emergence of an Efi- 
mov state of three magnons may be observed with ESR 
through an additional peak structure in the linewidth 
caused by the three-magnon resonance. Below we set 
?i = /cb = 1 and the unspecified range of integration is 
assumed to be from — oo to oo. 



II. QUANTUM CLUSTER EXPANSION 

Before we work on ESR in a dilute magnon gas, it is in- 
structive to demonstrate the quantum cluster expansion 
in a dilute Bose gas, which is described by a Hamiltonian 



H 



dr iplis^ 



iV 



A)Vr- 



dr^pl'lpl.'lpri'r- (1) 



Here Sp = /{2m) is the single-particle dispersion rela- 
tion, /i = — A < is the chemical potential, and the bare 
coupling constant g is related to the s-wave scattering 
length as by \/g — mA/(27r^) — m/(47ras) with A being 
a momentum cutoff. Our purpose in this section is to 
compute the single-particle spectral function 



A{uj,p) = -2ImG(a; -fiO+,p), 



(2) 
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where G{iujn,p) with a;„ = 27rn//3 being the Matsubara 
frequency is the Fourier transform of the imaginary-time 
propagator 

G{tu,,,p) = -iy^ I rfre'-"--^P- {Mr)^im- (3) 

As we will show in the next section, A{uj, 0) serves as the 
universal formula for ESR at low temperature and can be 
used to extract the scattering length between magnons. 
The imaginary-time propagator is defined by 

-iMrWO)) = -^Tv[e-P"Mr)^im, (4) 

where Z is the grand canonical partition function: Z = 
TT:[e~^^]. The systematic expansion over the fugacity 
z = e~^^ 1 can be developed by writing the grand 
canonical trace as a sum over canonical traces with fixed 
particle number N: Tr[-] = X]w=o ^^'^I ' ]■ Because 
trjv[e~'^"'^] oc , we obtain 



(r*(0)KzT) 



Z = l + V-^ + Oiz^) 



(5) 



where V is the system volume and A = ^J^iF^Jm, is the 
thermal de Broglie wavelength. Accordingly, the particle 
number density is found to be 



1 91nZ _ z 



+ 0(z2). 



(6) 



Therefore, the quantum cluster expansion is valid when 
the system is so dilute that the mean interparticle dis- 
tance is much larger than the thermal de Broglie wave- 
length: n"^/3 > A. 

The numerator in Eq. Q can be expanded over z in 
the same way. By denoting the Fourier transform of each 
term as 



y.tvM[e-^"Mr)4>lm-0{z'') 
the leading term is easily evaluated as 

1 — p-/3ep 



Go(ia;„,p) 



-Tz 

-Ep-A' 



(7) 



(8) 



On the other hand, after a straightforward calculation, 
the next-to-leading term is evaluated as 



Gl{iuJn:P) 

dq 



V/\^ -t- e 



F{iuj„,p; q) 



(2^)3 {luj,, 



Ay 



where 

F{iujn,p; q) = 



z + 0{z^), (9) 



iiT /m 



(10) 
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FIG. 1. Real and imaginary parts of a*{0)/{zT) in Eq. (fTi) 
as functions of A/os, represented by solid and dashed curves 
respectively. 



is the forward scattering amplitude between a particle 
with energy- momentum (ia;„ — A,p) and an on-shell 
particle with momentum q. Accordingly, the Fourier 
transform of the imaginary-time propagator ([3]) can be 
brought into the standard form 



G{iuJn,p) 



(11) 



where the quantum cluster expansion of the self-energy 
is found to be 



T,{iujn,p) 



dq 

(2^ 



■F{iujn,p;q)e-^''' +Oiz^). (12) 



Therefore, the self-energy at 0{z) is determined only by 
the two-particle physics, i.e, binary collisions with ther- 
mally excited particles. 

The resulting single-particle spectral function ^ at 
0(z°) is simply a delta function located a.t lu = Sp + A 
corresponding to the single-particle energy. By including 
the self-energy correction Y,{iujn,p) ~ 0{z) ^ 1, it be- 
comes a sharp peak whose line shape within the accuracy 
up to 0{z) is described by the Lorentzian 



^poak(i:^,P) 



-2Imcr(p) 



- Ep - A - Re (t(p)]2 + [Im a{p)]^ ' 

(13) 



where we introduced the on-shell self-energy 



a{p) = S(£p + A 
dq 



iO+,p) 
Stt/to 



~/3e<j 



(27r)3 l/as 



i\p- 



(14) 



Therefore, the energy shift and the decay width of a par- 
ticle with momentum p in a dilute Bose gas are given by 
the real and imaginary parts of a*{p), respectively. The 
corresponding normalized quantity a*{p)/{zT) ~ 0{1) 
at p = is plotted in Fig.[T]as a function of X/as- We find 
that the energy shift changes its sign and the decay width 
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FIG. 2. Single-particle spectral function A{ijJ,0) for z = 0.1 
and A/fls = 4 as a function of Suj = lj — A. The single- 
particle peak goes up to TA{ijj, 0) ~ 51 and there is a thresh- 
old singularity at the two-particle binding energy E2/T — 
-{\/as)'/{2n). 



reaches its maximum across the scattering resonance at 
ttg — )■ cxD. Such characteristic behaviors are universal and 
shared by ESR at low temperature. 

Besides the single-particle peak at w — Sp + A + 
Kea{p), there exists an additional structure in A{u!,p) 
a,t u! < Ep + A + E2 when two particles form a bound 
state with binding energy E2 — — 1/(7710^) < 0. This 
structure is due to the pole of the two-particle scattering 
amplitude (|T0)) and shows a square-root threshold singu- 
larity at uj = Ep + A + E2 in the single-particle spectral 
function as A{uj,p) ~ Z\/£p + A + E2 ~ uj. This feature 
at p = is demonstrated in Fig. [5] by choosing z — 0.1 
and A/fls = 4 as an example. 



system volume and the intensity of the incident radiation 
is given bjJ^ii^ 



■Imx(w + iO+,0), 



(16) 



where x(«w„,p) is the Fourier transform of the 
imaginary-time susceptibility 



XiiUJr 



Jo 



tpr 



{s~{T)s+m. (17) 



While the ESR experiments can measure the spectrum 
only at zero momentum, we develop the formulation for 
general p and set p = at the end. 

The ground state for a sufficiently large magnetic field 
B < is a fully polarized state with all spins point- 
ing downwards: S^\0) = -S\0) and S-\0) = 0. Ac- 
cordingly, we redefine the Hamiltonian to absorb the 
ground state energy Eo = -J2r(^JzS^ + DS^ - BS) 
so that H\0) = 0. Then at low temperature, magnons 
are thermally excited and thus the system becomes a di- 
lute magnon gas. A single magnon has the dispersion 
relation Sp — >S'J^-[1 — cos(p • e)] with the excitation 
energy A = -6SJ + 65 + 2SD - D - B. As long as 
the fugacity is small, z = e~^^ <C 1, the quantum clus- 
ter expansion can be developed for the dilute magnon gas 
similarly to the previous dilute Bose gas. The imaginary- 
time susceptibility is defined by 



1 



(18) 



where Z is the grand canonical partition function: Z = 
Tr[e~^^]. By writing the grand canonical trace as a sum 
over canonical traces with fixed magnon number N, we 
obtain 



III. ELECTRON SPIN RESONANCE 

A. Formulation 

We now apply the same method to ESR in a di- 
lute magnon gas. For definiteness, we consider a spin- 
s' Heisenberg model with exchange (J^ ^ J > 0) and 
single-ion {D ^ 0) anisotropics, which is described by 

H = —- (J S^S^^- + Jz S^S^^g) 



DY,is:r-Bj2s:. 



(15) 



Here J2e ^ J2e=±x ±y ±z ^ ^^'^ '^'^^^ vec- 
tors of a simple cubic lattice and spin operators = 
± iS^ and obey the usual commutation relations: 
[S+,S-] = 2S^Sr,r' and [S^,S$] = ±S^Sr,r'. ESR ex- 
periments measure an absorption intensity of electromag- 
netic radiation polarized perpendicular to the magnetic 
field axis. Within the linear response theory for a circular 
polarization, the absorption intensity normalized by the 



E 

N=0 



-I3H 



]^l + V- 



0{z\ 



(19) 



where we used trAr[e-^^] cx z and introduced the 
analogue of the thermal de Broglie wavelength by p = 
ae'^^^'^ /Io(2pSJ) with a being the lattice constant. Ac- 
cordingly, the magnetization relative to the ground state 
{Mq = —5") is found to be 



6M = 



1 d\nZ _ z_ 2 

y/3 dB " p3 + 



(20) 



The numerator in Eq. (jl8p can be expanded over z in 
the same way. By denoting the Fourier transform of each 
term as 



Xn 



{tuj„,p)= fdrY] e*-"^-*P- 

xtrjv[e"^^5-(T)5+(0)]^O(z^), 
the leading term is easily evaluated as 

1 - e-'^^pz 



Xo(«w„,p) = -2S- 



(21) 



(22) 
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On the other hand, after a straightforward calculation, 
the next-to-leading term is evaluated as 

Xl[lU}n,P) = -25 ^ T- 2: 

IW„ - Ep - A 

25 r' \/^.,/(^""'^''^L e-^--z + 0(z^), 



(23) 



where 



r(zc.„,p; g) = E [•^cos(H±2 . e) - J, cos(H^ • e)] 7(e) 

e 

- 2L>7(0) (24) 

is the forward scattering amplitude between a magnon 
with energy-momentum (ia;„ — A,p) and an on-shell 
magnon with momentum q. Here the unknown func- 
tion 7(r) = 7(— r) implicitly depends on {iuJn,P,Q) and 
satisfies the Lippmann-Schwinger equation 

7(r) = 2 cos(2^ ■r)+ . , , cos(fe ■ r) 

J-TT/al^Tr/a)'^ 



— tt/ a 

[^cos(£±2 • e) - J, cos(fe • e)] 7(e) - 20^(0) 



lUJn 



A + Eg - £{p+q)/2+k - £{p+q)/2-k 



(25) 



By setting r = x, y, z and r = 0, we obtain four cou- 
pled equations to determine 7(e) and 7(0) appearing in 
Eq. 

Accordingly, the Fourier transform of the imaginary- 
time susceptibility (jl7p can be brought into the standard 
form 



X(*w„,p) 



-2S 



iuJn - £p - A - ^{iu}n,p) ' 



(26) 



where the quantum cluster expansion of the self-energy 
is found to be 



^(iw„,p) = z 



1 1 a 



dq 



-tt/q (27r/a)' 



rr(zw„,p; q)e-^'-' +0(z^). 

(27) 



Therefore, the self-energy at 0(z) is determined only by 
the two-magnon physics, i.e, binary collisions with ther- 
mally excited magnons. The resulting ESR spectrum 
at 0{z'^) is simply a delta function located at a; = A 
corresponding to the single-magnon energy at p = 0. By 
including the self-energy correction S(ia;„,0) ^ 0{z) <C 
1, it becomes a sharp peak whose line shape within the 
accuracy up to 0(z) is described by the Lorentzian 



^peak('^) ~ -Z 



-2S'ImC(0) 



2 - A - Re^(0)]2 + [ImC(O)] 



(28) 



where we introduced the on-shell self-energy: C(p) = 
S{ep + A+iO~^ , p) . Therefore, the frequency shift and the 



linewidth of the single-magnon peak are given by the real 
and imaginary parts of ^*(0), respectively. Also, when 
two magnons form a bound state with binding energy 
E2 < 0, the ESR spectrum shows an additional struc- 
ture at w < A -I- i?2 similarly to Fig. [2l while we will not 
investigate it further. 



B. Solution 

Our remaining task is to solve the Lippmann- 
Schwinger equation (1251) with iujn = A + iO+ and p = 0: 



7(r) = 2 cos(| • r) -f- 



7T / a 



dk 



■ cos(fc • r) 



[Jcos(f ■ e) - J, cos(fc ■ e)] 7(e) - 2^7(0) 

£q - £q/2+k - £q/2~k + «0+ 



(29) 



to determine 7(e) and 7(0) in the two-magnon scattering 
amplitude (1^^ . In the low-temperature limit T ^ 0, the 
integration over q in the self-energy ((27|) is dominated by 
the region q ~ because of the Boltzmann factor. In 
this case, by expanding the right hand side of Eq. ((29)) 
up to 0(q), we find 7(0;) — ■y{y) — j{z) and the in- 
tegration over k can be performed analytically. Then, 
by substituting the obtained analytical solutions for 7(e) 
and 7(0) into Eq. ([24]) . we find that the two-magnon scat- 
tering amplitude takes the same form as that of bosons 
in Eq. ([T0| : 



lim r A + ?:0+,O;q = — — ' .. ,, 

q^o ^ ' l/as + i\q\/2 

= ^i?(A + zO+,0,q), (30) 

where l/m = 2SJa^ is the inverse effective mass of 
magnons and 

_ 2tt I-^ 3.7 .7 65,7 

a 25 - 1 + i(l - ^) + 3I^[1 ^)] 

(31) 

is the scattering length between magnons with 

dk 2 
(27r/a)3^,[l-cos(fc-e)] 

967r3 ^24/ V24/ V24/ V24/ ^ ' 

being the Watson's triple integral for a simple cubic 
lattice. We note that the same scattering length was 
obtained in Ref. 6 with a different approach. 

Accordingly, the low-temperature limit of the on-shell 
self-energy at p = appearing in Eq. (p8)) reduces to 
that in Eq. (ITil) . 



W : 



lim ^(0) 

T-)-0 



(7(0), 



(33a) 
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FIG. 3. Real (upper panels) and imaginary (lower panels) parts of ^*(0)/(2:r) in Eq. (|28|) for S — 1/2 (left panels), 5 = 1 
with D = (middle panels), and S = 1 with Jz = J (right panels) at T/J = 0.01, 0.1, 1, 10, represented by solid, dashed, 
dash-dotted, dotted curves, respectively. Two-magnon resonances as oo are located at Jz/J = 2.94, Jz/J — 4.87, and 
D/J — 4.77, respectively. 



and thus the hne shape of the single-magnon peak is de- 
scribed by the universal formula 

_lun /poakl^^) = ^S'^poak(w,0), (33b) 

where Apeak is the single- particle spectral function 
of bosons obtained in the previous section as Eq. (IT^ . 
This result is actually expected because the magnon gas 
at low temperature is so dilute that the system becomes 
independent of microscopic details and thus described by 
only a few low-energy parameters such as m and agJ^ 
Therefore, it is possible to extract the scattering length 
between magnons by fitting the universal formula (I33p to 
the experimentally measured temperature dependence of 
the line shape of the single-magnon peak. 

Away from the low-temperature limit, the line shape 
((28| is model-dependent and has to be evaluated nu- 
merically by solving Eq. (1^^ for general q. The fre- 
quency shift and the linewidth of the single-magnon peak 
are given by the real and imaginary parts of ^*(0), re- 
spectively, and the corresponding normalized quantity 
C{0)/{zT) - 0(1) is plotted in Fig. 13 For demon- 
stration, we choose three distinct cases where S — 1/2, 
S = 1 with D = as functions of Jz/J, and S = 1 with 
Jz = J as a function of D/J at four different temper- 
atures, T/J = 0.01, 0.1, 1, 10. Equation ^ indicates 
that the two-magnon resonances — >■ cx3 are located at 
Jz/J = 2.94, Jz/J = 4.87, and D/J = 4.77, respectively, 
where magnons interact strongly. We find that the line 
shape of the single-magnon peak is well described by the 
universal formula (j33l) at low temperature T < J and 
thus the frequency shift changes its sign and the linewidth 



reaches its maximum across the two-magnon resonance. 
Such characteristic behaviors become sharper with de- 
creasing temperature and can be seen moderately even 
at intermediate temperature T ~ J, while they disappear 
at higher temperature T > J. 



IV. CONCLUSION AND DISCUSSION 

In this article, we studied ESR in a dilute magnon 
gas that is realized in a ferromagnetic spin system at 
low temperature. We developed the quantum cluster 
expansion up to 0{z) which is determined by the two- 
magnon physics and showed that the frequency shift 
of the single-magnon peak changes its sign and the 
linewidth reaches its maximum across a scattering res- 
onance between magnons. Such characteristic behaviors 
are universal and can be used to experimentally locate 
the two-magnon resonance when an external parameter 
such as pressure is varied. Future achievement of the two- 
magnon resonance may have an impact comparable to the 
Feshbach resonance in ultracold atoms and will open up 
a rich variety of strongly correlated physics such as the 
recently proposed Efimov effect in quantum magnetsi^ 

It is straightforward in principle to continue this sys- 
tematic expansion to include higher-order corrections. In 
particular, the 0{z'^) term involves the three-magnon 
physics and thus it is possible to probe the Efimov ef- 
fect with ESR. When the system comes across the crit- 
ical coupling where an Efimov state of three magnons 
emerges from the scattering threshold, the linewidth of 



6 



the single-magnon peak as a function of the external pa- 
rameter is expected to show an additional peak struc- 
ture caused by the three-magnon resonance on either 
side of the two-magnon resonance. This feature is in 
analogy with ultracold atom experiments where an atom 
loss peak caused by the three-atom or atom-dimer reso- 
nance has been used as a signature of the emergence of 
an Efimov trimer^ii^ Similarly, a pair of universal four- 
magnon states associated with every Efimov statei^ii^ 
may be observed with ESR through an additional peak 
structure in the linewidth at 0{z^) caused by the four- 
magnon resonance. ^^'^^ Besides such characteristic be- 
haviors in the single-magnon peak, we expect an addi- 
tional structure in the ESR spectrum at cj < A -t- En 
when N magnons form a bound state with binding energy 
En < 0. This structure shows a threshold singularity at 
Lj = A + En as/(a;) - z^~^{A + En -cjY-^^-^I/^ , which 
may be used to measure binding energies of magnon Efi- 
mov states. Therefore, ESR is a powerful experimental 
technique to investigate the interaction among magnons 
and their spectrum. 

So far we considered ferromagnetic spin systems where 
the single-magnon dispersion relation has a minimum at 



zero momentum. On the other hand, it is also possi- 
ble to induce the two-magnon resonance and thus the 
magnon Efimov effect in spin systems with antiferromag- 
netic or frustrated exchange couplings where the single- 
magnon dispersion relation has a minimum at nonzero 
momentumi^i2£i2i In these cases, however, the sharp sig- 
nature of magnon scattering resonances discussed in this 
article will not appear in ESR because it measures the 
spectrum only at zero momentum. Therefore, different 
experimental techniques such as inelastic neutron scat- 
tering that can scan momentum space should be used to 
probe magnon scattering resonances. The quantum clus- 
ter expansion developed here will be useful to compute 
any relevant physical observables in a dilute magnon gas. 
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